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Coulomb screening and collective excitations in a graphene bilayer 

Xue-Feng Wang and Tapash Chakraborty 
Department of Physics and Astronomy, The University of Manitoba, Winnipeg, Canada, R3T 2N2 

We have investigated the Coulomb screening properties and collective excitations in a graphene 
bilayer. The static screening effect is anisotropic and is much stronger in the undoped graphene 
bilayer than in a monolayer graphene The dynamic screening shows the properties of a Dirac 
gas in the low frequency and that of a Fermi gas in the high frequency. The transition from the 
Dirac to the Fermi gas is also observed in the plasmon spectrum. Finally, we find that an electron 
gas in a doped graphene bilayer has quite similar properties as those of a Fermi gas in materials 
containing two energy valleys. 
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PACS numbers: 71.10.-w,75.10.Lp,75.70.Ak,71.70.Gm 



The realization of the monolayer and the bilayer 
graphene - the two-dimensional crystal of one and two 
layers of carbon atoms, has opened a new door for explo- 
ration of the fundamental physics and also fabrication of 
nanoelectronic devices [HlH. Being different from the 
multilayer graphene or graphite, the monolayer and the 
bilayer graphene are intrinsically zero-gap semiconduc- 
tors. The monolayer graphene has the properties of a 
chiral Dirac gas while the graphene bilayer has the energy 
band of a chiral Fermi gas at high energies with several 
quasi-Dirac points at the bottom of the band 0, [E IE 0] • 
Consequently, a comparison of the physical properties 
between them would offer new understanding and pro- 
vide interesting predictions about the different behaviors 
between the charged chiral-Dirac and the Fermi gases. 

In normal two-dimensional semiconductors, the many- 
body effects on the Coulomb interaction has been exten- 
sively studied using different methods. One of the most 
successful and widely accepted approaches is the random- 
phase approximation (RPA). In this approximation it is 
assumed that only the single-particle excitations of the 
same wavevector as the Coulomb interaction plays an ef- 
fective role in the screening process while the effects of 
others having different wavevectors cancel out. In the 
past years, there have been a few studies reported in the 
literature on the Coulomb screening and the collective ex- 
citation spectrum in the monolayer graphene using the 
RPA [1, [tj EE EL E2] • I n this paper, we present some of 
the interesting properties of the Coulomb screening and 
the collective excitations in a graphene bilayer. 

The graphene bilayer is formed by stacking two 
graphene layers in the same way as the stacking occurs 
in graphite, i.e., the Bernal stacking 0, H, 0, 0]. Each 
graphene layer has a hexagonal (honeycomb) carbon lat- 
tice [Fig.QJa)] which is composed of two periodic sublat- 
tices, A and B [Fig. QJb)]. In other words, there are two 
inequivalent lattice sites with atoms A and B in each unit 
cell of the periodic lattice. The two sublattices are dis- 
placed from each other along an edge of the hexagons by 
a distance of an = 1.42 A. In a graphene bilayer, there are 
four inequivalent sites in each unit cell, with atoms A and 
B at the top and A 1 and B' at the bottom. In the case of 
the Bernal stacking, the two graphene layers are arranged 



in such a way that the A sublattice is exactly on top of 
the sublattice B' with a vertical separation of 6o = 4 A 
[d| as shown in Fig. [Ha) and (b). The system can be 
described by the tight-binding model |5| characterized by 
three coupling parameters, 70 = 3.16 eV between atoms 
A and B or A' and B' (intra-layer coupling), 71 = 0.39 
eV between A and B' (the direct inter-layer coupling), 
and 73 = 0.315 eV between A' and B, A and A', or B 
and B' (the indirect inter-layer coupling). 

In the k space, the graphene bilayer has the same 
hexagonal Brillouin zone as that of a monolayer 
graphene. Its physical properties are mainly determined 
by the energy spectrum and the wavefunction near the 
two inequivalent corners of the Brillouin zone K and K' ', 
where the it* conduction band and 7r valence band meet 
at the Fermi surface [l4j]- Due to the strong interlayer 
coupling (the 7r orbit overlap) both the conduction band 
and the valence band in a bilayer are split by an energy 
of - 0.4 eV near the K and K' valleys 0(11 Gil. Since 
this energy splitting is larger than the energy range we 
are interested in from the bottom of the energy band, 
we take into account only the upper valence band and 
the lower conduction band. The bilayer graphene cannot 
be treated as two independent monolayer graphenes with 
the interlayer coupling as a perturbation because of the 
strong interlayer overlap of the it orbits. In contrast, the 
perturbation treatment of the interlayer coupling is valid 
for a normal double quantum well system [161 ] or in an 
intercalated graphite [r7| . 

In the effective-mass approximation, the electrons in 
the K valley is described by a Hamiltonian with a mix- 
ture of the linear and the quadratic terms of k 0, 0] 
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(1) 



with k± = k x ± k v and k = (k x , k y ) being measured from 
the K point. The effective mass of the quadratic term is 
= 2fi 2 7i/(3a 7o) 2 ~ 0.033 and the 'light' velocity of 
hk /2m* = 3a j 3 /2h w 10 5 m/s 
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the linear term is vq 
with fc » 10 8 /V3m _1 . 
The eigenfunction of 

, , , ih r ( p^fc 



the above 
with the 



Hamiltonian is 



energy 
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\h 2 k^fk 2 — 2fcofccos3<p + /2m* and the pseudospin 
angle \<f>k- Here <p — arg(fc + ), <p = &rg(ke~ 2tv — k^e 1 "-?) 
with arg(z) being the argument of a complex z = \z\e e , 
and A = 1(— 1) for the conduction (valence) band. 

For k ^> ko the electron states are chiral with = —2<p 
and have an approximately isotropic parabolic energy 
dispersion Efc = \h 2 k 2 /2m* . Near k — ko, the en- 
ergy dispersion becomes highly anisotropic as shown in 
Fig- [He) and (d) . The corresponding characteristic en- 
ergy is E a = kl/2m* = 3.9 meV. At E = 0, where 
the Fermi energy is located in undoped graphene bilayer, 
there are four contact points between the conduction and 
valence bands: One at k = 0, the center of the valley and 
three satellites at k = k in the directions of (p — 0, 7r/3, 
and 27r/3. They can be treated as four quasi-Dirac points 
because the electronic states near each point has a lin- 
ear energy dispersion and has the same chirality as those 
near a Dirac point in the monolayer graphene. However, 
compared to the monolayer grahene, the energy disper- 
sion here is anisotropic and the 'light' velocity is about 
ten times slower. As illustrated in Fig. []Jc) and (d), 
there is an energy pocket with a depth of about Eq/4 
at each quasi-Dirac point. The above peculiar charac- 
teristics of the graphene bilayer makes it quite different 
from the monolayer graphene in their Coulomb screening 
properties as described in the following. 

Just as for the case of the monolayer graphene [l(| and 
other SU(2) spin systems [lq ]. we find that the dielectric 
matrix of a graphene bilayer is again a unit matrix mul- 
tiplied by a dielectric function 

e(q,uj) = l-v q il (q,uj) (2) 

with the bare Coulomb interaction v q — e 2 /(2eoq) and 
the electron-hole propagator 

(q " j) + (3) 

The factor four comes from the degenerate two spins 
and two valleys at K and K 1 , f(x) is the Fermi func- 
tion, and the vertex factor reads \g^ X (<z)| 2 = [1 + 
AA' cos(#fc — Qk+q)]/2. Near the central quasi-Dirac point 
at k = 0, the intraband backward scattering and in- 
tcrband vertical Coulomb scattering are forbidden and 
| 5 ^- A (0)| 2 = \g^ x (~2k)\ 2 = 0. The same rules also hold 
for the three satellite quasi-Dirac points. For a large k 
(k > ko), |.9fe'~ A (0)| 2 = but \g x ' x (-2k)\ 2 = 1, i.e. the 
intraband backward transition is allowed but both the 
interband backward and vertical transitions are forbid- 
den. The above selection rules together with the energy 
dispersion of the carriers in the graphene bilayer indicate 
that electrons (holes) close to the bottom (top) of the 
conduction (valence) band have very different behaviors 
from those away from the bottom (top). 

Note that in a monolayer graphene the dielectric func- 
tion £ is invariant if all the parameters with the energy 
unit, to, Ep, and fcsT, and with the wavevector unit, k 




FIG. 1: (a) The hexagonal lattice of the top graphene layer 
(solid line) and the bottom layer (dotted line), (b) Periodic 
sublattices A and B' (dash-dotted lines), B (solid lines), and 
A 1 (dotted). In (a) and (b), the atoms of the sublattice A (B'), 
B, and A' are denoted by semi-filled, filled, and empty circles 
respectively, (c) The contour lines of the energy in the k x — k y 
plane near the K point. The corresponding energies, starting 
with the innermost curve are 0.1-Eo to Eq with an increment 
of 0.1i?o- (d) The energy spectrum for equally separated if 
from to 7r/6 (curves with increasing energy). 



and q, vary proportionally because of the linear energy 
dispersion of the Dirac gas [l(| • As a result, the dielec- 
tric function and the plasmon dispersion is uniform for 
systems with proportional parameters. In a graphene bi- 
layer, however, this is not true anymore because of the 
nonlinear energy dispersion. 

Coulomb screening: The static dielectric function at 
zero temperature versus q is plotted in Fig.^a). Its long 
wavelength limit is given by the properties of the four 
quasi-Dirac points. The central point has an isotropic 
'light' velocity vo = hko/(2m*) = 10 5 m/s while the 
satellite ones have the elliptic form of equienergy lines 
with a minimum 'light' velocity equal to vq along their 
radical direction and a maximum of 3vq along the az- 
imuthal direction. The static dielectric constant at q = 
is estimated to be e s = 1 + 3e 2 /(8eo^o) ~ 105. This 
value is much bigg er than the one for the monolayer 
graphene (4.5) [10| . This means that the long-range 
Coulomb interaction is much more strongly screened for 
the bilayer system, thanks to a much bigger density of 
states near the Fermi energy in a graphene bilayer. 

Another characteristic of the bilayer graphene is its 
screening anisotropy, especially for scattering at a dis- 
tance range of about 10 nm. This is shown by the differ- 
ence between the solid and the dotted curves in Fig. [21 
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FIG. 2: (a) The static dielectric function e a vs the wavevector 
q along the direction a — (solid) and tt/6 (dotted) at T — 0. 
(b) The same as (a) but at T = 4.2 K. (c) The real part of the 
dielectric function e r vs frequency ui at T — (dotted) and 
at T = 4.2 K (solid), (d) The imaginary part of the dielectric 
function Si vs ui at T = (dotted) and at T = 4.2 K (solid). 
In (c) and (d), q = 0.005 x 10 8 m" 1 and a = 0. In the limit 
— » oo, e r gradually approaches to one while e% approaches 
to zero. 
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FIG. 3: The plasmon spectrum of an undoped graphene bi- 
layer at T = (a) and at T = 4.2 K (b). The thick curves 
indicate the weakly Landau damped modes while the thin 
curves represent the strongly damped modes. 



corresponding to the directions of q pointing to any satel- 
lite quasi-Dirac points (a — 0) or to the middle of any 
two satellites (a = tt/6) respectively. Here a is the an- 
gle between q and the x-axis. At q = v3&o = 10 8 m _1 , 
the wavevector distance between any two satellite quasi- 
Dirac points, the anisotropy of e s reaches its maximum 
whit a mismatch of 20% along the different directions. 
The shoulder near q — ko — 0.58 x 10 s m _1 in the solid 
curve reflects the strong scattering between the carriers in 
the central and the tp — setellite quasi-Dirac points. At 
a finite temperature, the energy pockets near the quasi- 
Dirac points are partially occupied and the intraband 
scattering strength is greatly enhanced. As a result, the 
static dielectric function near q = increases rapidly, as 
shown in Fig. HJTd) at T = 4.2 K. 

In order to have a complete picture of the Coulomb 
screening and an insight for understanding the plasmon 
spectrum in bilayer graphene, we have calculated the real 
and imaginary parts of the dynamic dielectric function 
versus the energy for a small wavevector q = 0.005 along 
a = at T = and at T = 4.2 K in Fig. [2] (c) and (d). 

For oj > Eq/2, the dielectric function of a bilayer 
graphene is similar to that of a normal Fermi gas and 
its temperature dependence is weak. The step of e r and 
the peak of e% near ui = Eq/2 = 2 meV correspond to the 
single-particle excitation coupling states with a vanishing 



group velocity and having wavevectors located near the 
middle between the central and the satellite quasi-Dirac 
points. For small u>, however, the dielectric function be- 
comes more sensitive to the temperature and shows char- 
acteristics of the Dirac gas. One sign of the Dirac gas is 
the lack of Coulomb screening (e r w 1) in the energy 
window between 1 and 2 meV. Another sign is that a 
low-energy plasmon mode appears only at a finite tem- 
perature. As shown in Fig. [He), the e r has no negative 
value for the energy uj < Eo/2 at T — but evolves into 
a deep negative dip at a finite temperature T = 4.2 K, 
when the energy pockets near the quasi-Dirac points is 
partially occupied. As a result, one observes a weakly 
Landau damped plasmon mode of dispersion u> ~ ^fq at 
T = and a couple at finite temperatures. 

Collective excitation: In Fig. [3Ja), we calculate the 
plasmon spectrum of an intrinsic bilayer graphene (Ep = 
0) . The dispersion of the weakly Landau damped mode 
is indicated by the thick curve and has a ^fq dependence. 
Interestingly, the plasmon mode exists only at the energy 
higher than Eq/2, i.e., double the depth of the energy 
pockets in the quasi-Dirac points. At a finite temperature 
T = 4.2 K, another weakly damped plasmon mode shows 
up at the energy lower than Eq/2 and also has a disper- 
sion of ^/q near q = 0, as illustrated in Fig. Efb). The 
plasmon mode of higher energy existing at T = is not 
sensitive to the temperature. This temperature depen- 
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FIG. 4: The plasmon spectrum of a doped bilayer graphene 
(solid curve) with a typical carrier density of 10 12 cm" 2 . Cor- 
respondingly, E F = 36.3 meV and k F = 1.77 x 10 8 m _1 . The 
plasmon spectrum in the same system but without chirality is 
plotted as a dotted curve for comparison. Intra- (dark shaded) 
and inter- (light shaded) band single-particle continuums are 
also shown. 



dence of the low and high energy plasmon spectra repre- 
sents the marked difference between electron gases hav- 
ing linear (without collective excitation) and quadratic 
(with collective excitation) energy dispersion at T = 0. 
We have established earlier that the collective excitations 
appear only at a finite temperature for the Dirac gas [ic| . 
The electronic states in a graphene bilayer are similar to 
the Fermi type at the high energies but reverts to a Dirac 
type at the low energy. 

The carrier density of the system can be changed by 
doping 0. For a typical doping density of 10 12 cm~ 2 
[1,13' the Fermi energy is high enough from the conduc- 



tion bottom and the linear k term in the Hamiltonian can 
be neglected. The electrons then have a quadratic dis- 
persion but with chirality and <f> = —2ip. Near q = 0, the 
plasmon dispersion in a doped bilayer graphene has again 
a y/q dispersion as shown by the solid curve in Fig. [5] and 
shares the same dispersion uj 2 p D = [n e e 2 q/2eom*] 1 / 2 with 
a normal two-dimensional Fermi gas. To see the effect 
of the chirality, we plot as a dotted curve the plasmon 
dispersion of a normal two-dimensional Fermi gas with 
two valleys, for comparison. The two curves overlap for 
the small q but separate as q increases. The maximum 
difference in dispersion appears near q — ^/2kp when 
k and k + q form a right angle in the Fermi plane and 
the corresponding transition is forbidden in the graphene 
bilayer due to its chirality. 

In summary, we have explored the Coulomb screen- 
ing properties and the collective excitation modes in a 
graphene bilayer. In an undoped system, the static di- 
electric constant is much bigger than that in a monolayer 
graphene, due to the existence of four quasi-Dirac points 
in each energy valley and the much lower 'light' velocity 
of the quasi-Dirac points than the one in a monolayer 
graphene. The Coulomb screening also shows a strong 
anisotropy, especially near the wavevector q = \/3fco- 
The dynamic Coulomb screening has the characteristics 
of a Dirac gas on the low energy side and those of a 
Fermi gas on the high energy side. This transition from 
the Dirac to a Fermi gas is also reflected in the plasmon 
dispersion as the wavevector or the energy increases from 
zero. In a doped bilayer system, the long wavevector limit 
of the collective excitation has the same property as in a 
normal two-dimensional Fermi gas with two valleys. 
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